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Abstract 

The dynamics of a spherically symmetric thin shell with arbitrary rest 
mass and surface tension interacting with a central black hole is studied. A 
careful investigation of all classical solutions reveals that the value of the 
radius of the shell and of the radial velocity as an initial datum does not de¬ 
termine the motion of the shell; another configuration space must, therefore, 
be found. A different problem is that the shell Hamiltonians used in literature 
are complicated functions of momenta (non-local) and they are gauge depen¬ 
dent. To solve these problems, the existence is proved of a gauge invariant 
super-Hamiltonian that is quadratic in momenta and that generates the shell 
equations of motion. The true Hamiltonians are shown to follow from the 
super-Hamiltonian by a reduction procedure including a choice of gauge and 
solution of constraint; one important step in the proof is a lemma stating that 
the true Hamiltonians are uniquely determined (up to a canonical transfor¬ 
mation) by the equations of motion of the shell, the value of the total energy 
of the system, and the choice of time coordinate along the shell. As an exam¬ 
ple, the Kraus-Wilczek Hamiltonian is rederived from the super-Hamiltonian. 
The super-Hamiltonian coincides with that of a fictitious particle moving in a 
fixed two-dimensional Kruskal spacetime under the influence of two effective 
potentials. The pair consisting of a point of this spacetime and a unit time¬ 
like vector at the point, considered as an initial datum, determines a unique 
motion of the shell. 


1 Introduction 


One of the most spectacular and at the same time the most mysterious phenomena in 
nature is the gravitational collapse. It produces the largest fireworks in the sky and it 
ends up in a singularity, to which no known physical laws seem to be applicable. The 
enigma is enhanced by the appearance of horizons before the final stages of collapse. 
The horizon has some regularizing effect; for example, it prevents the divergence 
of the total energy due to concentration of matter (unlike the electrodynamics). 
However, the horizon leaks and there is another mystery: does the Hawking effect 
lead to violation of unitarity or not? 

The theory of gravitational collapse is difficult, not least because the phenomenon 
seems to belong to the high energy regime of quantum gravity and the progress in 
quantum gravity is slow. Thus, impatient people try to capture the core of the 
problem by working with simplihed models. This is also the line of the present 
paper. 

Let us consider the simplest models of the collapse: spherically symmetric thin 
shells (in general with a surface tension) falling in a held of a spherical black hole, and 
Oppenheimer-Snyder stars. Let us restrict the dynamical systems to their minimum: 
a hxed “amount of matter” without any internal degrees of freedom. For example, 
the total rest mass of dust shells and of the Oppenheimer-Snyder star is hxed. Let 
us, hnally, suppress all helds diherent from gravity. The resulting dynamical systems 
have just one degree of freedom; this can, at least locally, be described by the radius 
r of the matter part of the system (of the shell, of the Oppenheimer-Snyder star). 

Such or similar simple systems have been often utilized in literature. For instance, 
in the treatments investigating relativistic dynamics of spherical domain walls (some¬ 
times called bubbles) [|I|, models of gravitational collapse p, back reaction in the 
Hawking ehect [Q, or models of quantum black holes [Q, 0, 0. 

If we study these papers, we hnd rather surprisingly that the theory of even 
such simple systems is affiicted with severe problems. Thus, most papers start with 
guessing the Hamiltonian (or super-Hamiltonian) of the system to be studied from its 
equations of motion (eg. 0, i). However, it is well-known that a Hamiltonian 

and, therefore, an action is not uniquely determined by the equations of motion that 
it is to generate (this is the so-called inverse problem of variational calculus). In some 
papers, the Hamiltonians are derived from the hrst principles but the derivation is 
extremely complicated 0. The results of both guesses and derivations are, as a 
rule, very complicated, non-local Hamiltonians. A complicated Hamiltonian is a 
problem, if we are going to base a quantization on it and to calculate anything of 
interest. Sometimes, variables are chosen such that the corresponding canonical 
formalism becomes manageable but this simplicity is often bought by disregarding 
large portions of the space of solutions (the physical phase space); this seems to be 
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the case in |p. Finally, and this is really bad: different choices of time coordinate 
along the matter boundary lead to very different forms of Hamiltonian for the same 
system, thus making the corresponding naive quantizations non-equivalent, and so 
gauge dependent. 

The present paper is an attempt to deal with all these problems in the case of the 
shells. First, we show that a Hamiltonian of a system with one degree of freedom 
is uniquely determined (up to a canonical transformation) if its value is prescribed 
to be the total energy of the system. Very simple methods to calculate such Hamil¬ 
tonians are presented. Second, we present a careful study of all classical solutions 
of each respective system. This study reveals that the systems possess more states 
(ie. classical solutions) and admit more symmetry than various quantum approaches 
usually assume. There are not only bound states, but also scattering ones, and there 
is a number of states beyond horizon. The symmetries are related to isometries of 
Kruskal spacetime. This implies that the motion of the system is not (globally) 
determined by the value of the pair (r, r), and that the energy of the system is not 
a (globally) well-defined function of r and r. Another result is a one-to-one corre¬ 
spondence between the set of classical solutions and that of certain trajectories on a 
certain spacetime.Q The spacetime is closely associated with the classical solution; 
for the shell in the held of a black hole, it is the maximal extension of the uu-surface 
of the the part of Kruskal spacetime inside the shell and the trajectories coincide 
with those of the shell in this spacetime. The dynamics of the trajectories can be 
obtained from a super-Hamiltonian on which three conditions are imposed: it should 
be gauge invariant, at most quadratic in momenta and it should reproduce all sym¬ 
metries in the space of solutions. Under these conditions, the super-Hamiltonian is 
unique (up to an overal factor). In this way, the variables describing the gravita¬ 
tional field—which are pure gauges in the spherically symmetric case—are reduced 
away without any choice of coordinate conditions. A similar super-Hamiltonian has 
been obtained and used for the quantization of thin dust shells with the flat space- 
time inside |^. It follows from the uniqueness of the true Hamiltonians that they 
can be obtained from our super-Hamiltonian by reduction procedures. 

The plan of the paper is as follows. In Sec. ^ we briefly describe the relevant 
properties of Kruskal spacetime, in particular the symmetries. In Sec. |], we prove the 
uniqueness of Hamiltonian discussed above, present some simple methods to derive 
such Hamiltonians and show two examples of complicated and gauge-dependent 
Hamiltonians. Sec. ^ is devoted to massive thin shells with arbitrary surface tension. 
By including tension, we can describe shells made from fundamental fields such as 

^S. K. Blau, E. I. Guendelman and A. H. Guth following the suggestion of S. Coleman, 
discussed the classical solution of the equation of motion of the spherical domain wall as the 
solution of the equation of motion of a particle moving in one dimension under the influence of a 
potential. 
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scalar fields in the bubbles arising in cosmological phase transitions Sec. ^ 
collects what we need about null shells. In Sec. |^, we write down a Hamiltonian 
formalism for the shells. In Sec. starting from our super-Hamiltonian, we derive 
the Hamiltonians that describe the motion of null shells in the Kruskal spacetime 
external to the shell (like Ref. |^). 

2 Kruskal spacetime 

All generic (not extreme) black hole spacetimes contain a region which has the same 
topology as the Kruskal spacetime. It includes two identical stationary asymptot¬ 
ically flat submanifolds that are causally separated by two crossing horizons, and 
two identical dynamical submanifolds that are time inversions of each other. For 
the astrophysical description of a stellar collapse, only a half of this spacetime is 
necessary and relevant. What can be the consequence of the existence of the other 
half? Apart from the discussions of primordial wormholes or white holes, it seems 
that it may have some importance in the quantum theory: indeed, for example the 
popular explanation of Hawking effect uses the states beyond the horizon because 
they have negative energy with respect to inhnity. Let us start by briefly describing 
the Kruskal spacetime, restricting ourselves to the two-dimensional version of it. 

Consider with double null coordinates u and v and the metric 

ds^ = — —^e~^{—dudv), ( 1 ) 

F 

where i? > 0, G is the Newton constant and F = F{—uv) is dehned by its inverse, 

F~^{x) := {x — l)e'^ 

in the interval x G (0, cxd) and —uv G (—l,oo). The metric is well-dehned and 
analytic in the region —uv G (—1, cxo). This part of together with the metric (|ID 
is called Kruskal spacetime (M,g). 

The Kruskal spacetime is time and space orientable. We will call future the 
orientation in which the function [u + v) is increasing, and right in which [v — 
u) is increasing. There are two independent discrete isometries that invert these 
orientations: T, dehned by m i—>■ —n, v i—>■ —u and C by m i—>■ n, n i—^ m. Their 
composition I = CT = TC is called the Kruskal inversion] I inverts both orientations. 
Observe that C is a new transformation, enabled by the topology of the Kruskal 
manifold. C has nothing to do with parity; indeed, the parity transformation P 
leaves both u and v unchanged but transforms d and ip hy d h-^ tt/ 2 — d and 
ip ^ ip + TT in the full four-dimensional manifold. 
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There is also a continuous one-dimensional group of isometries generated by the 
Killing vector field 


e: = 


AGE 


d d ^ 


If g{s) is an arbitrary element of this group, then 


( 2 ) 


T^(s)T = ^(-s), Cg{s)C = g{-s), lg{s)l = g{s). 


The signature and orientation of ^ changes over M. Let us define four quadrants in 
M as follows: 


Qi 

:= {(m,-!;) 

M < 0, n > 0}, 

Qii 

:= {(m,-!;) 

\u > 0,v < 0}, 

Qiii 

:= {(m,-!;) 

\u > 0,v > 0}, 

Qiv 

:= {(M,n) 

M < 0, n < 0}. 


Then, the orientation and signature of ^ is given by the table 



signature 

orientation 

Qi 

timelike 

future 

Qii 

timelike 

past 

Qiii 

spacelike 

right 

Qiv 

spacelike 

left 


The group generated by T, C, and ^ is called the Kruskal group, Gk- 
Two important functions, t and r, on M are defined as follows: 


r := 2GE E{—uv) 


everywhere. 


t := AGE arctanh 


u + V 
—u + V 


in Qj U Qjj, and 

r-. 1 

t := AGE arctanh- 

u + V 

in Qni U Qiv- The pair (t,r) of functions can be considered as coordinates within 
each of the four quandrants; these are usual Schwarzschild coordinates. We have 
the following transformations: 


on M, and 


-.-F 


dr = —2GE——{ydu + udv) 
F 


dt = -2GE 


vdu — udv 


uv 


A 




on Qj U Qii, and 


vdu — udv 


dt = 2GE 

on Qijj U Qiv, where we can also set 


uv 


F{—uv) 


r 

2GE’ 


—uv = 


r-2GE 

2GE 


r 

q2GE . 


Within the quadrants and in the Schwarzschild coordinates, we have ^ = d/dt. 

The behaviour of the functions t and r with respect to the Kruskal group is as 
follows. Let p & M and g{s) be an element of the component of identity of Gk- 
Then 


{t{g{s)p),r{g{s)p)) 

(t(Tp),r(Tp)) 

{t{Gp),r{Gp)) 

{t(lp),r(lp)) 


{t{p) + s,r{p)), 
{-t{p),r{p)), 
{-t{p),r{p)), 
{t{p),r{p)). 


Thus, r is an invariant function of Gk, and any invariant function of Gk must be a 
function of r. Moreover, on the part of M covered by the four quadrants, the hxed 
value (t,r) of the pair of functions t and r dehnes exactly two points, p and I(p). 


3 Hamiltonians of simple systems 

In this section, we consider spherically symmetric systems that a) are sufficiently 
simple so that they have just one degree of freedom and b) the matter has a well 
dehned external boundary so that the geometry outside this boundary coincides with 
a portion of the Kruskal spacetime. We study the ways, in which a Hamiltonian 
generating the motion of such systems can be obtained from hrst principles (that 
is, from the Einstein-Hilbert action to which some matter action is coupled). Our 
aim is to show that such Hamiltonians are a) unique in certain sense if the time 
coordinate is chosen, b) very complicated, and c) gauge dependent. 

In order to obtain a true Hamiltonian, we have to reduce the system. This, 
roughly, is the following procedure. First, we have to adapt the coordinates to 
the spherical symmetry. Second, we have to choose a foliation given by levels of a 
function T on the phase space and a radial coordinate X along the folii; then, we 
have to solve the constraints for the momenta pt and px conjugate to T and X. 
The result can be written as 


Pt + H{q,p,T) = 0, 
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where q and p are the coordinates on the physical phase space of the system; the 
true Hamiltonian is H{q,p,T) (see, eg. ||T0|). 

Since the spacetime outside the matter has the Kruskal geometry, we can choose 
the foliation such that it assymptotically coincides with the Schwarzschild foliation. 
It is well-known (0, 0) that the value of the true Hamiltonian then coincides 
with the ADM energy E [E = GM, where M is the Schwarzschild mass parameter 
of the external Kruskal spacetime). Observe that the foliation can be different from 
the foliation by the Schwarzschild time coordinate at the boundary of the matter. 
Thus, the foliation will determine a time parametrization of the boundary. Let us 
denote the time parameter by s; the value s{p) of the parameter at a given point 
p of the boundary coincides with the value T{p) of the function T at p. In fact, s 
can be an arbitrary time parameter along the boundary. In this way, each motion 
of the system determines a function r(s), where r is the Schwarzschild radius of 
the boundary at the value s of the time parameter. As the ADM energy is always 
conserved, H cannot depend explicitly on T. 

Since the system has just one degree of freedom, the function r(s) can contain 
all information about the motion (at least locally). In this case, we can choose r 
and its conjugate momentum p as the canonical coordinates on an open set in the 
physical phase space of the system. The reduction procedure will then result in the 
true Hamiltonian of the form El (r, p). 

On the other hand, one can also study the dynamical equations of the system. 
They follow from Einstein’s equations and from the stress-energy tensor conservation 
for the matter. Solving the gravity part of these equations by the Schwarzschild 
metric, a non-trivial radial equation remains: 


r = p{r,E), 


(3) 


where the dot denotes the derivative with respect to the parameter s. The system 
can possess different sectors, so p can also depend on a few discrete parameters. 
We stress that the function p is completely determined by the dynamical equations; 
to derive its form, one does not need to employ the canonical formalism and its 
reductions. 

The crucial question now is: How many Hamiltonians of the form H{r,p) will 
lead to the equation of motion (j^)? It is well-known that there are many such 
Hamiltonians . However, if we impose the condition on H that its value coincides 
with the total energy E, then the possible Hamiltonians are strongly limited. Indeed, 
suppose that H (r, p) is such a Hamiltonian. Then, we must have 

. dH 
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This equation represents the transformation between the variables (r, r) and {r,p). 
Thus, the Eq. implies that 

This can be considered as a differential equation for hf; it is even a separable hrst- 
order ordinary differential equation as r plays the role of a parameter: 


dp = 


dH 


p{r,Hy 


( 4 ) 


Hence, the general solution has the form: 


P = Po{r) + J 


dH 


p{r,Hy 


( 5 ) 


where po{r) is an arbitrary function of r. Solving Eq. (^) for H, we obtain a family 
of Hamiltonians of the form: 


H = H{r,p - po{r)). 

The result is clearly non-unique, but this non-uniqueness, represented by the func¬ 
tion ^ 0(^)5 can be regarded as the symmetry of the formalism under the transfor¬ 
mation 

r I— >■ r, p p — pf){r). (6) 

It is clear that we always have the freedom to perform such a canonical transforma¬ 
tion in the classical Hamiltonian formalism (to simplify the equations, for example). 
More important, however, is the fact that the transformation (|^) belongs to the small 
class of canonical transformations that can be implemented by a unitary transfor¬ 
mation in the quantum theory. Indeed, suppose that we represent the Heisenberg 
commutation relations for r and p as usual: 

(fiV^i)(r) = rV’i(r), (pi^/>i)(r) = 

and that the wave functions etc. span the Hilbert space with the scalar 

product 

poo 

('0i,(/)i)=/ dry)l{r)(t)i{r). 

JO 

Next, we consider another copy of the same Hilbert space, but denote its elements 
by etc. and dehne a unitary transformation between the spaces by 

y)2{r) = 
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where Y (r) is a fixed differentiable function. The form of the operators r 2 and p 2 is 
given by 

f2 = = fi, 

p2 = = Pi — y' ■ 

Hence, the choice Y{r) = f drpo{r) shows our claim. In this sense, the Hamiltonian 
H can be considered as uniquely determined by 1) the equations of motion, 2) the 
value of energy and 3) the gauge. 

The description of dynamics based on a true Hamiltonian has, however, the two 
disadvantages discussed in the Introduction: complicated form and gauge depen¬ 
dence. 

Let us consider a null dust shell (see Sec. |); it always moves along a null geodesic 
in both Kruskal spacetimes, the internal as well as the external one with respect to 
the shell. In the Schwarzschild coordinates, Eq. ([T0|) holds, and so the functions t+(s) 
and r(s) that describe the trajectory of the shell in the external Kruskal spacetime 
must satisfy the equation 

-f+{r)t+^ + fp\r)r^ = 0 . 


Hence, 


and 


' dr 
dt^ 


p(r, E) = e 


= 1 - 


2GE\‘ 


1 - 


2GE 


where e = ±1 distinguishes the ingoing from the outgoing shells. In this case, Eq. 
(^ becomes 


r 


P = Po- In 


1 - 


2GE 


and its solution with respect to E yields the form of the Hamiltonian: 

^-e—(p-po(r)) 
r 




(7) 


This is a very complicated, non-local Hamiltonian. 

Next, we choose the Kraus-Wilczek coordinates r and r on the external Kruskal 
spacetime which are defined by H 


v/f- 72^ 


V2^ 


r = f_|_ -|- 2V2GEr + 2GE In 
so that the metric takes the spatially fiat form 

= —fj^dr"^ -|- 2\l ^ drdr + dr'^ -|- r'^dGP. 


( 8 ) 


r 



















The equation of a radial null geodesic in this metric reads 


drV l2GE dr 2GE 

Tt) + ~ 


0 


and its solution is 


dr l2GE 



where e has the same meaning as above. The Eq. now reads 


Gp = Gpo{r) - V2^ r — er In 


\/2GE — e^/r 

a/t 


(9) 


This will coincide with the Kraus-Wilczek momentum p^ if one matches the units 
and the notation and if Po(^) is chosen properly: 

Gpo{r) = sj2GE_r + er In 


V2^EZ-e^ 

Jr 


where E_ is the Schwarzschild energy of the Kruskal spacetime inside the shell. We 
observe again that the Hamiltonian is very complicated (there is no ecxplicit formula 
this time that would express the solution of the Eq. with respect to E), and it 
is also very different from Eq. 0 . 

In order to solve this problem, we shall try to hnd a super-Hamiltonian, which will 
be gauge independent, at most quadratic in momenta, and such that its reduction 
based on a suitable gauge will give each of the complicated Hamiltonians. Such a 
super-Hamiltonian can either be itself directly used for quantization, or dehne the 
factor ordering of the complicated Hamiltonians. 


4 Massive shells 

In this section, a careful study of the states of massive spherically symmetric dust 
shells (in general with arbitrary surface tension) will be given. In the following 
section, the null dust shells will emerge as a (degenerate) simple limit case. 

4.1 The trajectory of the shell 

Consider two Kruskal manifolds Ati and M .2 with energies E_ and E+ containing 
isometric timelike hyper-surfaces Si and S 2 which divide each Atj in two parts, 
left M.i-, and right Afi+. Let us solder Afi- to M. 2 + along Si and S 2 . Most 
points in the resulting surface S := Si fl S 2 in the resulting spacetime will have a 
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neighbourhood in which we can introduce Schwarzschild coordinates on both sides. 
The metric left (—) and the metric right (+) is given by (e = ±1) 


dsl = —f^{r)dt‘^ + ^{r)dr‘^ + r‘^{d9‘^ + sin^ 9 dcp'^), ( 10 ) 

where 


we assume that r, 9 and 0 are continuous functions through S, however, in general, 
at S. The surface S is given by 

fe = fe(s), r = r(s), 9 = id, 0 = 9 ?, (11) 

where s, ■d and ip are coordinates on E and s is the proper time along the radial 

generators of S. Let us denote by 

Tt ■= (4,r, 0,0) 

the tangent vector along these generators. Since s is the proper time, we obtain 

fete = T,^fe + r‘^, ( 12 ) 

where := sign(/£0). The meaning of the left hand side is fJe = where 

is the Killing vector (|^). This would be a total energy with respect to the right 

inhnity for a particle moving in the respective Kruskal spacetime, and it would be 
a conserved quantity for a geodesic motion. Let be the normal vector to E, 
being external to M. 2 +, and m_ to Afi-; we assume to be normalized, or 

fdem° - = 0 , 

fe{rrfif - f:\rn\f =-1. 

It follows that 

m\ = + (13) 

Je 

where := sign(m(‘0^r). 

There is an important relation between a and r that we are going to derive (we 
omit all epsilons in the indices). The Killing vector is tangential to r = const 
curves; hence, is orthogonal to the vectors dehned by := Thus, 

there are positive numbers ^ and p, almost everywhere along E, such that 





(the timelike vector first) is an orthonormal dyad. Similarly, is an or¬ 

thonormal dyad at the same point so that 

= a^ + b^, 

^ P 

Km^ = b— + a —, 

^ P 

where a and b are some reals and k, = ±1; if k = 1), the two dyads have the 

same (opposite) orientation. Then, 

= -a(sign^^)^, 

cr = sign{g^i,m>^p'') = sign(Ka(signp'")), 
where signn^ = -|-1 (—1), if is spacelike (timelike). Thus, 

a = K(signp^)(sign^^)(-r) = kt. 


because (signp^)(sign^^) = —1. k can easily be determined. We distinguish two 
cases: is future or past oriented, and introduce the corresponding sign rj = ±1. 

Consider hrst the future oriented case. In Qj, and are both timelike and 
future oriented, so a > 0. and are both spacelike, p^ is oriented to the right 
and em^ to the left; hence, = —e. In Qji, p^ and are both timelike and both 
future oriented, so 6 > 0. and are spacelike, oriented right and em^ left; 
again, = —e. A similar analysis for the remaining quadrants conhrms the relation 
(Te = —CTe for future oriented shells. Suppose now that is past oriented. Then 
r just changes sign, but everything else (in particular, the a’s) remains as before; 
hence, (Te = er^ for past oriented shells. Thus, the results can be summarized by the 
equation 

cr^ = -pen. (14) 

Finally, we calculate the acceleration of the shell in the spacetime directed 
out of M^. It is dehned by 

a, = m^pV^TP, 

where 

XI'^TP •= TP T rpflrpl/ 

^ s-^- e • -^e ' ^ e e e ") 

and TP^^ are the Christoffel symbols of the Kruskal metric gefiu- Substituting for A 
from Eq. (|T^) and using Eq. (|T^, we obtain after some rearrangements 


a, = a. 


(//: 


(15) 


These relations will be important for the exact counting of possible solutions of 
Israel’s equation, which determines the motion of the shell. 
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4.2 The matter of the shell 


We assume the shell to be made of ideal fluid. A spherically symmetric shell matter 
is “isotropic”, and so it is always “ideal fluid”, though the fluid could have some 
internal degrees of freedom (cf. |^), which we are suppresing. Thus, we set 


= (p + p)T^r + p7 


kl 


where 7 ^^ := is the induced metric on S, 

_ dxi^ 

ek- 

y"^) the embedding functions of S in Atip. If we choose functions (|TID 
as embedding functions we obtain 

7 ^^ = - 1 , 7 ^^ = r^(s), = r^(s) sin^ 


p is the surface mass density and p is the surface pressure (negative of surface 
tension). We assume that the state of the shell matter depends only on its total 
surface area A := An A and that the kind of the material is dehned by the barotropic 
equation of state p = p{p). We will also assume that the stress-energy in the 
shell is conserved (which follows, if there is vacuum outside the shell, from Israel’s 
equation—see later). Thus, 

T^\i = 0 , 

where the vertical bar denotes the covariant derivative associated with the metric 
7 fc/. If one substitutes for the stress-energy tensor, one obtains the energy equation 


(pT%+pT% = 0. 


This is equivalent to 


+ p + p(fj) = 0 . 


(16) 


The solution to this equation will be of the form p{A) and it will depend on one 
constant. We dehne our system to be just one particular solution to Eq. (^). We 
assume that the solution dehnes a positive density p for all r G (0, cx)). Thus, 
there is a well-behaved unique rest mass density at any radius. In this way, the 
system is as elementary as possible: it has no internal degree of freedom and its 
total energy depends only on its radius and velocity. Finally, we dehne the so-called 
mass function M{r) by 

M(r) := A{r)p{A{r)); 


this dehnition is justihed by the fact that M(r) will appear in many equations at a 
prominent place (in some cases, it represents the proper mass of the shell). 
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As an example, consider the eqnation of state 


P = kp, 

where k G (—1, +1). The correponding solntion to Eq. ([l^) reads 

p = 


where mi is a constant. If fc > 0 (positive pressnre), then M(0) = oo and M{oo) = 0; 
if A; < 0 (positive snrface tension), then M(0) =0 and M{oo) = cx); for fc = 0 (dnst), 
M = const. 


We snmmarize: 
form 


the non-zero components of the stress-energy tensor have the 


Tss ^ 


'I?!? 





P{P) 

sin^ P 


(17) 


4.3 Israel’s equation 


The dynamics of a (massive) shell is governed by equations found by Dautcourt 
and brought to a nice geometrical form by Israel |K|. Israel’s equation can be 


written in the form 


Qkl ^ Qkl ^ 

where Qf := — Lf and L^ki '■= ~'kkie^i-,vX^kXei is the second fundamental form 

of S in Ail- or At 2 +- 

The way (|l^) of writing the Israel equation makes the independence of the equa¬ 
tion on the orientation of the surface S manifest. It is also invariant with respect 
to any coordinate change of y's, in particular with respect to the time orientation 
change i—>■ —y^. If we use the embedding formulas (|T2), we obtain 


L aa 




= fe + r'^, 


= -ro-eV/,-f-r^sin ^ 


where is given by Eq. (lA). Then the trace reads 


and 


Qr 


gr 


L. = —a. — 2a, 


Vfe + 


= 2a, 


Vfe + 


1 f , Vfe + r^' 

— r a, -F (Je - 

rj^ ^ \ ly' 


sin^ P 


1 ( , V/e + r 

- ' Qe + ae - 


r 
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If we substitute these formulas into the L.H.S. and Eqs. (0 into the R.H.S. of 
Eq. (p!8[), we obtain Israel’s equation for our model: 



= 4:7iGrp, 

(19) 

\Jf+ + r^ 

Qj-\- Cl— “1“ (T_|_ ~h C7— 

= — SvrGp. 

(20) 


If we substitute for from Eq. (0 into then the resulting equation together 
with the time derivative of Eq. (|^ implies the energy equation (PIBD, or the time 
derivative of Eq. (W) together with Eq. (IR) implies the resulting equation. 

We can summarize: the full dynamics of the shell is contained in the following 
equations: 


fie = -r]eae\Jfe + r‘^, 

( 21 ) 


( 22 ) 


where the meaning of the paremeters r/, e, cXe (all equal to ±1) was explained above, 
/g = 1 — 2GEe/r, and M(r) is the mass function dehned in Sec. [4.2| . 


4.4 The radial equation 

Double squaring Eq. (|2^), we obtain the radial equation 

r2 + v{r) = 0, (23) 


where 


V{r) = 


G^M^r) G{E+ + E_) {E+- E_f 




+ 1 . 


(24) 


From the point of view of the existence and properties of solutions, Eq. (^^ is easier 
to study than Eq. (^2]). Let us do this study hrst and postpone the question of how 
the solutions of these two equations are exactly related. 

Eq. (p3|) is invariant with respect to the transformations 


s <->■ —s, or s H->• s + So, 

and 

E+ E_, or M —M. 

Suppose that the pair (E_,E+) is given. Eq. (|^) together with the initial con¬ 
dition r(0) = tq has either no solution, if E(ro) > 0, or two solutions, r(s) and 
r(—s), if E(ro) < 0. Thus, for each maximal interval (ri,r 2 ) in which V{r) < 0, we 
obtain just one or just two maximal solutions, depending on what happens at the 
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points ri and r2, where V{ri) = V{r 2 ) = 0. If V'{r 2 ) 7^ 0, then we can extend the 
solution r(s) by r(—s) at this point; similarly, if ri 7^ 0 and V'{ri) 7^ 0. In each of 
these cases, there is only one solution in the interval. If none of these cases takes 
place, then we have two maximal solutions, r{s) and r{—s), in the interval; they 
can be chosen in some standart way for each interval. Any other solution within 
the interval can be obtained by the shift s 1—^ s + Sq so that it satishes the initial 
condition r(so) = ’"o for any Tq G (ri,r2). The solutions that reach inhnity can only 
exist if {Ei,E 2 ) is such that Id(cx)) < 0. 

Let us consider an example: the dust matter, which has p = 0 and thus M = 
const. For the potential (|^) we obtain 


V{r) = 


4 



G(B+ + B_)(i) 


(g+ - E-f 
M2 


+ 1 . 


Then 


so that 




dV 


in the whole range r G (0, cx)). Further, 

1/(0) = -CX), I/(cx) = 1 


{E+ - E_) 
M2 


For given E^ and there is a unique such that V{rm) = 0 (allowing for the 
value Tm = C)o), and the range of possible initial values rg for which a solution exists 
is then [0,rm]- Hence, we always have a unique solution satisfying r(0) = 0, if 
Tm < 00. There are two solutions, r(s) and r-(—s), satisfying r(0) = 0 if = cxo. 

We shall have hyperbolic solutions if V( 00 ) < 0, a parabolic solution if V( 00 ) = 0, 
and elliptic solutions if H(cx)) > 0. For the elliptic (bound) solutions, < 00 . The 
unique non-negative solution of the equation V(r) = 0 is 

The parabolic solutions satisfy £'+ — E^ = ±M. If the maximal radius reaches a 
horizon, = 2GE^, the condition 


E., - E, 


M2 

4^ 


must be satished. 

We can easily understand the spectrum of E^ for a hxed E_. Three cases can be 
recognized according to the value of i?_: 
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Case < M/2 


If G (-E- + M, cx)), we have the hyperbolic solution s, if = E_ + M parabolic, 
and for Ej^ G (0,i?_ + M) elliptic. At Ej^ = {E_ + ^E'^ + M‘^)/2, = 2GE^; if 

E+ decreases beyond this value, then also decreases but it does not reach 2GE_; 
the minimal value of is (GM^)/(2(M — E_)). 

Case E_ G [M/2, M) 

This is the same as the previous case but if E^ decreases after reaches the value 
2GT^+, Tm can also reach the value 2GE _—at Ej^ = E_ — M'^/{4E_)] then, as Ej^ 
decreases further, increases again reaching the value (GM^)/(2(M — E_)) > 
2GE_ at E+ = 0. 

Case E_ G [M, cx)) 

The hyperbolic solutions fill two intervals: for E^ G (£'_ + M, cx) they reach the 
right inhnity, and for E^ G (0, E- — M) they reach the left one. We have two 
parabolic solutions, £'+ = E_ ± M, and in the interval E^ G {E_ — M, E_ + M) the 
solutions are elliptic. 

4.5 Relation between solutions of the radial equation and 
the shell spacetime 

The radial equation ( P3|) is a consequence of Israel’s equation (P^). Hence, every so¬ 
lution of Eq. (p^ must solve Eq. (p3[). However, as (^) is obtained by squaring (p2[) 
twice, some solutions to (p3D need not solve (|22|): one could loose the information 
about the signs cXg. 

For the understanding of this problem, the following observation is crucial. If one 
calculates the expression from the radial equation (^), one obtains: 



( 26 ) 


( 26 ) 


Thus, the expression on the left hand side of Israel’s equation. 



(27) 


with some signs a± equals to 


CT+ I J+1 -|- (T_ I J_ 
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It follows immediately that the only choice of the signs that satishes Israel’s equation 

(0)is 

(Tg = signJe. (28) 

To obtain a valid shell spacetime, we have to satisfy also Eqs. (0) and (0; this 
leads to 

feir{s)% = -T]eJ,. (29) 

Thus, quite surprisingly, each solution of the radial equation (p3|) generates a solution 
of Eqs. ( ^TD and (p^ ) which describe the full dynamics of the shell. 

Let us construct a shell spacetime with a future oriented shell motion from a 
solution of the radial equation: 

1. Choose an oriented pair {Ei, E 2 ) of positive numbers; this dehnes two Kruskal 
spacetimes, and we will construct a shell spacetime by having M.i with Ei = 
E- to the left and M .2 with E 2 = E_|_ to the right. 

2. Choose a solution r(s) to the radial equation. Then the expressions given 
by Eq. (pSD are well-dehned for each s and e. The function r(s) solves then 
Eq. (p^ with the signs a^. 

3. Eq. (p9|) determines the functions t^{s) up to an additive constant, te(s) h->■ 

te{s) + Q. 

4. The pair of functions r{s),t^{s) given on some common interval of s determines 
two different curves in each Kruskal spacetime A4i and A42- They are related 
by the inversion map I (cf. Sec. |^), so only one of them is future oriented. Thus, 
in each of A4i and M 2 , there is a unique surface Si and S2, respectively. We 
cut the spacetimes along it and solder the parts together. The geometry of the 
resulting spacetime does not depend on q: the constants can be made zero by 
continuous isometries in each half. 

The construction also revealed some ambiguities. If we did not £x the orien¬ 
tation of the shell movement and the order of the pair (E_,E^), we would have 
obtained four non isometrical spacetimes in general. This corresponds to the num¬ 
ber of elements in the discrete subgroup of the Kruskal group Gk- Reversing the 
order of (E_,E+) corresponds to the map C (cf. Sec. |^) performed in each space- 
time (All, El) and (M 2 , E 2 )—the dependence of r on s along any curve cannot be 
changed by such a map. Similarly, if we allowed also for past oriented shells, then 
the two functions t^(s) and r(s) would determine two curves in each spacetime that 
are related by the map I; by inverting the parameter s h-^ —s, we obtain a future- 
oriented shell motion. Let us denote the parameter inversion by T^. It is clear from 
this discussion that a description of the dynamics that is based solely on the radial 
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trajectory r(s) will not distinguish physically different states. More concretely: the 
radial equation p3| ) is quadratic in the total energy if considered as an equation 
for it has two roots in general. Hence, Ej^ is not a well-defined function on the 
space of allowed pairs (r, r). 

In fact, there are even more states to be distinguished: sometimes, two spacetimes 
that are otherwise completely isometric to each other can describe different physical 
states. The crucial point is whether the isometry moves also the observers or not. 
Let us assume that the observers in our model are assembled at the right spacelike 
inhnity. If we move the shell by the continuous subgroup of Gk, t+ ^-i- + c+, 
without moving the observers, then the physical state will change: the observer will 
be able to detect c+ as a change of the shell arrival time at a fixed radius, for 
example. Thus, one and the same radial trajectory r(s) represents 4 x cx) different 
physical states (or 2 x oo, if one does not like the past-oriented shells). Observe that 
the time shift t- t- + c_ does not do anything to change the physical state with 
respect to the right family of observers. Thus, we have to couple the internal and 
the external times, if we wish to describe the shell dynamics by means of variables 
that refer to the left (internal) spacetime. 

Another important point is that the validity of our equations is not limited to 
spacetimes with just one shell; as far as there are no crossings, there can be other 
shells and/or sources. Interesting examples of two-shell spacetimes can be easily 
constructed using the symmetries. Let us take three Kruskal spacetimes, (AIi,i?i), 
{A 42 , E 2 ), and (Ada, E 3 ) such that Ei = E 3 . The shell E12 between Adi and Ad2 can 
be taken as the C- or loTg-image of the shell S23 between Ad 2 and Ads. Consider a 
source lying far left in (Adi, Ei) so that both shells are to the right from it. Such a 
source can be constructed explicitly by taking for example a suitable Oppenheimer- 
Snyder star. The motion of both shells is not changed by this source. We can 
interpret this as follows. The difference ei2 := E 2 — Ei is the contribution of the 
energy of the shell S12 to the total energy of the system, and similarly 623 := E^ — E 2 
for S23; these contributions exactly cancel each other: 


ei2 — —623- 


One can consider e as the gravitational charge of the shell. Then S12 and S23 have 
opposite gravitational charges. It is interesting to observe that the opposite charges 
give rise to opposite accellerations. Let us calculate the acceleration ai_ of E12 in 
Adi directed out of Adi and the acceleration —03+ of S23 in Ads directed into Ads 
(both pointing right). Using the formulas (pAj), (pS]), (p^) and (^), we obtain an 
important formula for a^: 

-±T 
— , Je- 
dr 
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Substituting ei 2 or 623 , respectively, for in Eq. ( p 6 D we obtain: 

_ ^( GM{r) ei 2 \ 
dr y 2 r M{r) J 

^ d / GM(r) ^ 623 \ 
dr y 2 r M(r) j 

However, ei 2 = — 623 , so ai_ = —(— 03 +) as claimed. 

An interesting example is the following. Let S 23 be a hyperbolic shell outgoing 
to the right of the event horizons of both spacetimes (M 2 ,-^ 2 ) and (^^ 3 ,^ 3 ); it 
will reach the right assymptotic region sometimes in future. The other shell, S 12 , 
can be constructed by the map loT^ of the hrst one in (M 2 ,E 2 ) and (^^ 3 ,^ 3 ), 
where the map I is interpreted as a map I ■. M 3 ^ Mi. Thus, S 12 is a hyperbolic 
outgoing (with respect to the right inhnity) state lying left from the event horizons 
of both (Ml, El) and (M 2 ,E 2 ). If we could construct a quantum theory in such 
a way that these will be a pair of shells whose classical interpretation will coincide 
with our construction, and if we can get creation of such pairs, spontaneous or 
induced by some agent, then this may give a model of Hawking effect with some 
sort of back reaction (cf. P]). There seem to be other possibilities, too: S 12 can be 
constructed from S 23 by C. One should try to implement these speculations by a 
suitable quantization of the system. 


5 Null shells 


The spherically symmetric shell from null matter moving in vacuum has a very simple 
theory |^, []^. Its motion with respect to each of the embedding spacetimes Mi 


and M 2 is that of radial null geodesics. Thus, if a = 0,1, are two coordinates in 
the Kruskal un-plane, and geab(x) the corresponding metric, the dynamical equation 
of the shell takes the form 

= 0 , 

where the dot is the derivative with respect to an arbitrary parameter (in a two- 
dimensional spacetime, every smooth null curve is a null geodesic). In local Schwarz- 
schild coordinates (pT]) we obtain 


- uti + = 0 . 


(30) 


There are only six different types of null geodesics in any two-dimensional Kruskal 
spacetime (we consider only the future oriented shells): 

1 . generic geodesics, that reach the right (left) asymptotic region; we distinguish 
them by the sign parameter = -|-1 ((^ = — 1 ). 
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2 . outgoing (Cf > 0) and ingoing {(r < 0) for each C-tyP® (out- and ingoing is 
defined with respect to the right inhnity). 


3. outgoing and ingoing horizons. 


Along any generic null geodesics, the coordinate r is an affine parameter and it 
behaves strictly monotonically, acquiring all values from the interval (0, cxo). There 
are only two exceptions: the horizons. 

Given two generic null geodesics of the same type, one in Afi and the other in 
A42, they dehne the surfaces Si and S 2 , and there is only one soldering of Afi_ 
with M 2 + along Si and S 2 : the function r must be continuous. This is possible for 
any values of B_ and In particular, the outgoing (ingoing) horizon in Afi can 
be soldered only to the outgoing (ingoing) horizon in Ad 2 , and this only if Ei = E 2 . 
However, for a null dust shell (which is pressureless), the soldering must be affine 
(the shell is “affinely conciliable”—cf. 0), and this is not unique. For the outgoing 
horizon, using the double null Kruskal coordinates u and n, a general soldering is 
given by U- = 0, = 0, = n_ -f a, where a is an arbitrary constant, and 

similarly for the ingoing case: n_ = = 0, «+ = «_-!-/?,/? arbitrary. 

It follows that we can dehne the quantity ex for all affine parameters A along any 
null shell. We dehne hrst cae by 




where is the Killing vector (^ of the respective Kruskal spacetime. As x“(A) is a 
null geodesic in each of the Kruskal spacetimes, cae is conserved along the curve on 
both sides. If the shell reaches an asymptotic region, both geodesics reach the same 
(right or left) asymptotic region in the respective spacetime, and so both eA+ and 
ca- have the same sign. From Eq. (|30|) , we obtain 

f ■ 

^"dXj ~ \dx) ’ 



but, in the Schwarzschild coordinates, 

_ r dte 


(32) 


Hence, eA+ = ex-- Moreover, if A = r, then exe = C- For rmll shells, we thus have 
ex- = ex+, and we can dehne 


e\ ■= e\- = ex+. 

The structure of matter of the given null shell constructed according to the recipe 
above is given by the Barrabes-Israel equations. In particular, for a generic shell, 
Eqs. (52) and (53) of |]^ imply 
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1 ) the pressure p inside the shell vanishes (we have null dust if the shell is sur¬ 
rounded by vacuum), 

2 ) the shell stress-energy tensor is given by 


rpab _ /- laib 

- ^ 47rr2 ’ 


(33) 


where is the tangent vector to the radial null geodesics affinely parametrized by 
r; in the shell coordinates y^, 

dyk 


f = 


dr 


For a horizon shell (see [T^), the only non-vanishing components of the stress- 
energy tensor with respect to the coordinates {v,'d,ip) (or {u,'d,ip)) are 

e e 


rpVV _ 


-a, (T““ = 


:/ 5 ), 


(34) 


{2>2GE)H ’ ^ {2>2GE)H' 

where e is the basis of natural logarithms and E = E- = Ej^. 

The energy density of the shell is a measurable quantity (see m- In order that 
this be positive, the condition 


E+-E.= C|B+ - E. 


(35) 


must hold for generic shells, and some conditions on the sign of a and (3 for the 
horizon shells. 

From the formulas (^OD, (|33|) and (^41) it follows that the matter content of the 
null shell is independent of their trajectories; they seem to have an internal degree 
of freedom. However, it is possible, at least formally, to create a relation between r 
and Ej^ — E_ by hand, using the arbitrariness of an affine parameter A. This will 
turn out to be necessary in the Hamiltonian formalism where the hnal justihcation 
for this step lies. 

Let us dehne: the affine parameter A along shells is a physical parameter, if 

(36) 


where, as in the case of massive shells, r] = -|-1(—1) if E is future (past) oriented. 
(This also keeps the mass density positive for generic shells.) Then, Eqs. (pT]) and 
(0) imply 

This is the “radial equation for the null shells” corresponding to the Eqs. (PBf ) and 
(P^) for massive shells. Together with the time equation 



vC 


dr 

dX 


(38) 
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it determines the trajectory (the left hand side does not depend on e). 

In this form, the null-shell dynamics is just a limit case of the massive-shell 
dynamics (with the exception of the horizon shells). Indeed, let us dehne the physical 
parameter A for massive shells by 


(] s 


Then, the radial equation (^) and {‘M) becomes 


T) 


(39) 


where 


V{r) = 


G^M^ir) GM‘^{r){E+ + E_) 


{E+ - E_f + M^r). 


(40) 


Clearly, hmM(r)=o V{r) = —(T^+ — E_Y^ and Eq. (|39|) goes over into Eq. (^). 


6 Shells as particles on two-dimensional Kruskal 
spacetimes 


In this section, we will try to rewrite the dynamics of the shells as a dynamics of a 
hctitious particle on a fixed two-dimensional Kruskal background. 

Let us study the motion of a point particle on a two-dimensional Kruskal space- 
time with the energy E that results from the action 

S = J ds{pox^ + pix^ — AfH), (41) 


where are some coordinates on the spacetime, which serve as canonical coordi¬ 
nates of the particle, pa are the conjugate momenta. A/" is a Lagrange multiplier. 


n 


9‘^(Pa - UiMPi, - U(i) + 


(42) 


is the super-Hamiltonian, g°'^{x) is the Kruskal metric in the coordinates U{r[x)) 
and W{r{x)) are some functions. Eq. (|4^) provides the most general (up to an overall 
factor that is independent of the momenta) super-Hamiltonian Td with the properties 


1 . is a polynomial of second order in momenta, 

2. the momentum dependence of 7d has the correct light-cone structure, 

3. EC as a function of x^ and p^ is invariant with respect to Gk- 
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The action is clearly invariant with respect to all transformations of coordinates 

Within each quadrant of the Kruskal spacetime, we can transform the canonical 
coordinates to Schwarzschild coordinates t and r. The resulting action is 

S = J ds{pti + PrT — AfH), (43) 


where 


n 


-r\pt + uff + fpl + w^ 


(44) 


and /(r) := 1 — {2GE)/r. 

The action (^5D is meaningful only inside the four quadrants but can be considered 
as a result of a canonical transformation from the variables which can be global 
Kruskal coordinates u and v. The solutions can be matched through the horizons 
and the matching is well-defined by the tranformation to the Kruskal or Eddington- 
Finkelstein coordinates. 

Variation with respect to pt and pr yields: 


ft = -Af{pt + Uf), 
r = Affpr- 


We suppose hrst that W 0 and choose the gauge J\f = p/W, where the sign p 
corresponds, as before, to different time orientations of the shell; then 


Pt + Uf = -pWfi, (45) 

Pr = pWf~^r, (46) 

and the constraint (following from the variation of the action with respect to A/") 


n = 0 (47) 

implies that the parameter s is a proper time. Eqs. (0-13). and the equation that 
results from the variation of the action with respect to t, 


Pt = const. 


(48) 


form a complete system of equations of motion. The equation that is obtained by 
varying the action with respect to the variable r is a consequence of the system 

we obtain 


Inserting Eqs. (^6]) and (^8]) into Eq 


f + r^ = 


' Pt + Uf 
. W 


(49) 
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It is easy to see that Eqs. ( ^91) and (|i5| ) are equivalent to the radial equations (p5D, 
(p^) and the time equation (]29|), provided that the constants E and pt are chosen 
properly. Let E = E_. Then Eq. (PB|) with e = —1 will be satished if 


Pt 

= -E+, 

(60) 

Uf- 

^ GM\r) 

2 r ’ 

(61) 

w 

= M(r). 

(52) 


The apparent singularity in the super-Hamiltonian Ti can be removed if Ti is mul¬ 
tiplied by /(r) because provides another factor /(r) at 

There are two important elements in the choices 

1. The choice E = E _: this makes the shell to move on the left spacetime and thus 
decouples the background along which the shell moves from the gravitational 
held the shell produces. A great simplihcation is a consequence. 

2. Eq. (|50|) : we require that the total energy of the full system, which is E+, 
coincides with this conserved quantity. 


The above considerations enable us now to formulate the following theorem. 


Theorem 1 There is one-to-one correspondence between the states of the motion of 
the shell with mass function M{r) and the dynamical trajectories of a point particle 
with Pt < 0 which follow from the variational principle with the action 


S' = y ds{paX^ — MTC), 


where 


1 r 


^ = 2 [9-\Pa - U-{r)U{Pb - U-{r)^t) + M\r) 


g_{x) is the Kruskal metric with as energy parameter, and 


UJr) = 


f- 


E_ - 


GM^{r 

2 r 


(53) 

(54) 


(55) 


Proof The construction of the shell spacetime from a given solution of the 
variational principle (p3D goes as follows. The super-Hamiltonian is chosen such 
that the solution is identical with the shell trajectory in the spacetime (M,E) = 
(Ml, E_). Then, setting E+ = —pt, the spacetime (M 2 , E^) is well-dehned and we 
have to hnd the shell trajectory there. The function r(s) is the same as in Mi, but 
t+(s) must satisfy 


f+i+ 


/GM(r) 
V 2r 


E+-E \ 
M(r) J ■ 
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By this the trajectory in A ^2 is hxed up to the map given by the discrete isometry C 
and a map given by any element of the component of identity of the Kruskal group 
Gk (see Sec. The map by C does not change the resulting spacetime because 
always only one and the same half M. 2 + of M .2 can be soldered to M.i-. The time 
shift can be hxed by some coupling of the internal time to the external time 
For example, we can require that the point of the trajectory in M.i with = 0 
coincides with the point of the trajectory in M .2 with f_|_ = 0. (Another example is 
given in the next section.) We thus obtain a unique state of the shell. On the other 
hand, given a shell spacetime, the trajectory of the shell in (Adi, EJ) is well-dehned 
and it is a solution of the variational principle (|53|). There is indeed a one-to-one 
correspondence, Q.E.D. 

Let us show that the variational principle (^B] ) yields also the dynamics of the null 
shells. If we set M{r) =0 in Eqs. (|5^) and (^31), we obtain for the super-Hamiltonian 
in Schwarzschild coordinates 

^ = ^[-/T^(pt + E_)2 + /_p2]. (56) 

Let us choose the gauge Af = rj. Then variations with respect to pt and pr yield 

f-i- = -p{pt + E_), (57) 

r = pf-Pr- 

Inserting this into the constraint = 0, we obtain 

-f-il + fZ^r^ = 0. 

Thus, we have a null geodesic. Moreover, the comparison of Eqs. (^) with Eq. 
(p^) shows that s is the physical parameter now. The trajectory of the shell in 
{M. 2 i E+ = —pt) is determined similarly as for the massive case; the soldering, 
however, is unique only for the generic shells. For horizon shells, an affine shift 
stays undetermined; and only this shift determines the spacetime and the energy 
density of the shell (cf. Eq. (^). 

What is the use of the super-Hamiltonian (^)? The discussion of Sec. allows 
only the following accurate statement. Any true Hamiltonian H that results from 
Ti by some reduction procedure is a meaningful one, if its value coincides with the 
total Energy of the system (that is with E^). Those reduction procedures of the 
constraint system (^) that result in a different value of H are not allowed. 

The things have been arranged in such a way that the allowed reduction proce¬ 
dures are for example those which use a stationary foliation. That is, given any 
spacelike surface S of the auxiliary Kruskal spacetime, the stationary foliation is 
dehned as the family {g{s)S \ s G R}, where g{s) is an element of the component 
of identity of the Kruskal group Gk- By these foliations, however, not all allowed 
reduction procedures are exhausted, as an example in the next section will show. 
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7 Out-Side Story 


In this section, we shall transform the dynamics of the shell to the form of particle 
dynamics on the external spacetime {A42,E^). Snch schemes are in use but they 
are complicated. Our aim is to show that they can be obtained from the variational 
principle (|53[) by a reduction procedure including a choice of gauge. We shall restrict 
ourselves to the null shells. As an example, we shall rederive the Kraus-Wilczek 
Hamiltonian [Q from our super-Hamiltonian. 

The physical phase space (whose points are maximal classical solutions) of the 
generic null shell consists of four (topologically separated) sectors.^ The subsequent 
analysis applies to the sector of those shells that reach the right infinity and are 
outgoing. The other sectors can be dealt with in a completely analogous manner. 

To describe the dynamics of the outgoing shells, we choose the Eddington-Fin- 
kelstein retarded coordinates u and r. The metric in the internal spacetime reads 

ds^ = —f-du^ — 2dudr -|- r‘^dVL^. 

In these coordinates the super-Hamiltonian (Q) with M(r) = 0 is given by 

n = -{pu + E_){pr + fZ^E.) + i/_(p, + fZ^E_)\ 

but we allow only the solutions with p^ = —fZ^E_ as the others, with pu = {f-Pr — 
E_)/2, belong to the ingoing sectors. Thus, we can simplify the super-Hamiltonian 
to Pr + fZ^E-. The equations of motion then read 

ft = 0, r =J\f, Pu = ^, Pr = -fZ^E_. 


The phase functions u and pu are thus perennials (constants of motion) spanning 
the physical phase space f in this sector (cf. [0); on T, the perennials have the 
ranges 

uE{—oo,oo), PuE{—oo,0). 


The constraint surface T is the submanifold Pr = —fZ^E_ of the extended phase 
space T with the coordinates r, pr, u and pu- Thus, we can cover the constraint 
surface T by the coordinates r, u and pu {r E (0, cxo)). 

The reparametrization invariant system defined above is to be reduced so that a 
“true” Hamiltonian results. For that purpose, first a family of transversal surfaces 
{Ti}—the so-called time levels—is needed (see |T^ for more detail): T* C T for each 
t, and Tt have exactly one (transversal) intersection with each c-orbit (trajectory) 


^The reader is referred to Ref. ||^ and Ref. 0 for some reviews of the constrained dynamics 
of reparametrization invariant systems in which the concepts we use in this section are described 
in detail. 
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given by M = const, pu = const, r = A. (Here t is any suitable time parameter—not 
necessarily the Schwarzschild time). can be defined by an equation of the form 

t = T{r,u,pu), 


where T{x, y, z) is a suitable function of the three variables. is a symplectic space 
for each t] as a canonical chart on F^, {u\rt,Pu\rt) can be chosen. Second, for a 
construction of a Hamiltonian dynamics, a family of time shifts 6tt' '■ Fj/ i—^ Fj is 
necessary such that each 6tt' is a symplectic diffeomorphism |^]. Clearly, in our case 
a trivial family of 9tt' is obtained by mapping the points {u,Pu) € F^/ to {u,Pu) G F^; 
this would give a trivial (frozen) dynamics. A more suitable possibility is to employ 
the momentum pu as a generator of time shifts. Since 


{r,Pu} = ^, {m,P„} = 1, 

and (r, u) can be considered as Eddington-Finkelstein coordinates, Pu generates just 
the Schwarzschild time shift. The requirement that Fi _|_5 is obtained from Fj by Pu5 
leads to the following condition on T: 


{T,pJ} = 5 V5gR, 


or 


The general solution has the form 


dT 

du 


= 1 . 


T{r,u,pu) = u + T{r,pu). (58) 

The next requirement we impose is that the transversal surfaces can be interpreted 
as time-surfaces in the external (right) spacetime. To find such surfaces, we have 
to couple the external time to the internal one in such a way that the time in 
can be expressed as a function of our phase space variables. The most 
straightforward and useful way is to require that the retarded time be continuous 
across the shell (this can easily be achieved since it is determined up to a constant), 

M = = M+. (59) 


Then, using E_|_ = —pu, we can specify the function T. For example, if we want 
to employ the surfaces of constant Schwarzschild time then we use the formula 
expressing in terms of u+ and r. 


4 


M+ -|- r -|- A In 


r — A 
A 


(60) 
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where A := —2Gpu- Eq. (|^ ) then yields the following form of T (satisfying the 
condition (^8|)): 

T{r,u,pu) = u + r + A In — -j— . 

Finally, instead of the perennials u and pu, we would like to write our dynamics 
by means of the radial variable r and a corresponding conjugate momentum p. Let 
us study this transformation for general T. Of course, T must be such that the 
function r is a non-trivial coordinate along L*, that is, there must exist a solution 
r = R{t — u,pu) of the equation t = u + T{r,pu). Thus, we have the identity 

t-u = f{R{t - u,Pu),Pu)- (61) 

To hnd the momentum p{t,u,pu) canonically conjugate to the coordinate r = R{t — 
u,pu), we look for a canonical transformation (for each t) 


{u,pu) ^ (r,p). 

The transformation will be canonical, if 

dR dp dR dp 
du dpu dpu du 

Let us search for p in the form 


(62) 


p = P{R{t - U,Pu),Pu)- 


Then Eq. (^) is equivalent to 

dP \ 

dpu)^ 




From the identity (|UT|), we obtain the hnal result for the momentum conjugate to r 
in the form 

P{r,pu) = Po{r) - J (63) 

where Po{r) is an arbitrary function of r. 

The Hamiltonian of Kraus and Wilczek. Let us use the formula (|63|) to rederive 
the Kraus-Wilczek Hamiltonian obtained in 0. Kraus and Wilczek employ the 
system of coordinates r and r in the external Kruskal spacetime (Af 2 , E+) that has 
been dehned in Sec. ^ (formula (||)). The relations (|60| ) and (H) imply that the 
retarded time can be expressed by means of r as follows: 



u+ = T — r — 2vC4r — 2y4 In 











Thus, the transversal surfaces are dehned by the function T of the form: 


T{r,u,pu) 


u + r + 2a/]^ + 2^4 In 



1 


The function T satishes the condition (|58|) because the foliation of the external 
spacetime (Af 2 , -£-+) by the Kraus-Wilczek time r is stationary. Comparing the last 
relation with Eq. (^), we hnd T, from which we obtain 


df 

dr 


1 



Hence, the formula (|63D becomes identical with Eq. (^) with e = -|-1 and —pu = H. 

We can conclude: the non-locality of the Hamiltonian of Kraus and Wilczek is 
not due to any physical property of the system, but to the choice of gauge and to 
the reduction procedure. 
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